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I. INTRODUCTION 

Simultaneous heat and mass transfer from different geometries embedded in porous media has many 
engineering and geophysical application such as geothermal reservoirs, drying of porous solids, thermal 
insulation, enhanced oil recovery, packed - bed catalytic reactors, cooling of nuclear reactors, and underground 
energy transport. Convection in porous media has applications in geothermal energy recovery, oil extraction, 
thermal energy storage and flow through filtering devices. Convective heat transfer in porous media has 
received considerable attention in recent years owing to its importance in various technological applications 
such as fiber and granular insulation, electronic system cooling, cool combustors, and porous material 
regenerative heat exchangers. The phenomena of mass transfer are also very common in theory of stellar 
structure. Viscous mechanical dissipation effects are important in geophysical flows and also in certain 
industrial operations and are usually characterized by the Eckert number. 

Soundalgekar V.M. et al [1], was studied on flow past oscillating plate with variable temperature. 
Soundalgekar V.M. et al [2], was discussed free convection effects on the flow past a vertical oscillating plate. 
Vajravelu K., et al [3], presented a solution for Heat transfer in a viscous fluid over a stretching sheet with 
viscous dissipation and internal heat generation. Pop I, et al [4], was analyzed viscous dissipation effects on 
unsteady free convection flow past an infinite plate with constant suction and heat source. Tania S. Khaleque, et 
al [5], was studied Effects of radiation, heat generation and viscous dissipation on MHD free convection flow 
along a stretching sheet. Mansour M.A., et al [6], was investigated radiative and free convection effects on the 
oscillatory flow past a vertical plate. Soundalgekar V.M, et al [7], presented a solution for Effects of free 
convection currents and mass transfer on the flow past a vertical oscillating plate. Soundalgekar V.M, et al [8], 
discussed Effects of mass transfer on the flow past an oscillating infinite vertical plate with constant heat flux. 
Bala siddulu malga , et al [9], was analyzed viscous dissipation effect on unsteady free convection and mass 
transfer flow past an accelerated vertical porous plate with suction. Hitesh Kumar [10], was investigated 
radiative Heat Transfer with Hydro magnetic flow and viscous dissipation over a stretching surface in the 
presence of variable heat flux. Israel - cookey, C., Ogulu, A., Omubo - Pepple, V.M.,[11], was studied The 
influence of viscous dissipation and radiation on unsteady MHD free convection flow past an infinite heated 
vertical plate in a porous medium with time dependent suction. Mahajan, R.L., Gebhart, B.B.[12] presented a 
solution for viscous dissipation effects in Buoyancy- Induced flows. Zueco Jordan, J., [13], was discussed 
Network Simulation Method Applied to Radiation and Dissipation Effects on MHD Unsteady Free Convection 
over Vertical Porous Plate. Anjana Bhattacharya et al, [14] were analyzed Theoretical study of Chemical 
Reaction Effects on Vertical Oscillating Plate Immersed in a Stably Stratified Fluid. Saras wat Amit and 
Srivastava R.K.,[15] was discussed Heat and Mass Transfer Effects on Flow past an Oscillating Infinite vertical 
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plate with Variable Temperature through Porous Media. P.M. Kishore, et al,[16] was studied Effects of Heat 
transfer and viscous dissipation on MHD free convection flow past an exponentially accelerated vertical plate 
with variable temperature. V.Rajesh, [17] was investigated Heat Source and Mass Transfer Effects on MHD 
Flow an Elastic -Viscous Fluid through a porous medium. 

The objective of the present paper is to study the unsteady flow of an incompressible viscous fluid past 
an oscillating vertical plate on taking into account of viscous dissipation, heat and mass in presence of variable 
temperature. The dimensionless governing equations are solved by using an implicit finite difference method of 
Crank - Nicolson’s type. 



II. FORMULATION OF THE PROBLEM 

Consider the unsteady flow of an incompressible viscous fluid which is initially at rest past an infinite 
vertical plate with variable temperature through a porous medium. The flow is assumed to be in x-direction, 
which is taken along the vertical plate in the upward direction .The y-axis is taken be normal to the plate. 
Initially the plate and the fluid are at same Temperature with same concentration level at all points. At 
time t' > 0 the plate starts oscillating in its own plane with a velocity u = u 0 cos co't'. The plate Temperature is 
raised to T w and the levels of concentration near the plate are raised to C w linearly with time t. The viscosity is 
taken into account with constant permeability of porous medium. Then by usual Boussinesq’s approximation the 
unsteady flow is governed by the following equation. 




du 



d z u 



ir = gP( T' - T'oo ) + gp (C - C 00 )+ u -Ej - v(rr) 



dt 
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(1) 

( 2 ) 
( 3 ) 



with the following initial and boundary conditions 
t'<0, u = 0, T = T „o C’ = C, for ally, 
t' > 0, u = u 0 cos co't', T =T o0 + (T w -T 00 )At 
C = Coo + ( C w - Coo ) A t , at y = 0. 
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u = 0, T — > C — > C oo as y— >oo. t' > 0. (4) 

where A = ^ ^ 



On introducing the following dimensionless quantities 
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Initial and boundary conditions in non dimensional form are 



U = 0, 


0=0, 


C = 0 , 


for all 


Y < 0, 


t<0 




U = cos cot, 


0 =t, 


C = t, 


at 


II 
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(9) 



III. METHOD OF SOLUTION 

Numerical technique: Equations (6) - (8) are coupled non-linear partial differential equations and are to 
be solved under the initial and boundary conditions of equations (9). However exact or approximate solutions 
are not possible for this set of equations and hence we use implicit finite difference technique of Crank- 
Nicolson method. The finite difference equations corresponding to equations (6) - (8) are as follows 
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Initial and boundary conditions are in the following form 
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Here the suffix i corresponds to Y, j corresponds to t and L— ► oo. Also At = t j+1 — 



tj and 



AY = Y i+1 — Yj knowing the values of C,0,u, at a time t, we can calculate the values at a time t + At as 

follows we substitute i = 1, 2, 3, N-l in equation (12) which results in a tridiagonal system of 

equations in unknown values of C. using initial and boundary conditions the system can be solved by Thomas 
algorithm. Thus C is known at all values of Y at a time t + At. Then the knowing the values of C and applying 
the same procedure and using boundary conditions we calculate © from equation (11) the knowing the values of 
© and applying the same procedure and using boundary conditions we calculate u from equation (10). This 
procedure is continued to obtain the solution till desired time t. 



IV. RESULT AND DISCUSSION 

The velocity, Temperature and concentration profiles have been computed by using Crank- Nicolson 
finite difference method. The numerical calculations are carried out for the effect of the flow parameters such as 
Phase angle ( cot = tt/ 4 ), Thermal Grashof Number ( Gr = 2 ), Modified Grashof Number ( Gc = 5 ), Prandtl 
Number ( Pr = 0.71(air) Pr = 7(water) ), Schmidt Number ( Sc = 0.6 ), Eckert Number ( Ec = 0.5 ), 
Permeabilityparameter ( K = 0.5 ), time ( t = 0.2 ). The velocity, temperature and concentration profiles has 
been studied and presented in fig 2 to 9. 

❖ Effect of Phase angle (□'$): The velocity profile for different values of phase angle 

( cot = 0, 7i/6, 7i/4, 7i/2 ) is presented in Figure 2. It is observed from the figure that the velocity 

decreases with increasing phase angle ( cot ) 

❖ Effect of Thermal Grashof Number ( Gr ) : The velocity profile for different values of thermal Grashof 

Number ( Gr = - 10 , - 5 , 5 , 10 ) is presented in Figure 3. It is observed from the figure that the velocity 

increases with increasing Thermal Grashof Number ( Gr ). 



♦> Effect of Modified Grashof Number ( Gc ) : : The velocity profile for different values of modified 
Grashof Number (Gc = -5,-3, 3 , 5) is presented in Figure 4. It is observed from the figure that The 
velocity increases with increasing modified Grashof Number ( Gc ). 

❖ Effect of permeability parameter ( K ) : For the different values of permeability (K = 5, 10, 15 ) 

the velocity profile is drawn and is presented in Figure 5. It is observed from the figure that the velocity 
decreases with increasing permeability parameter ( K ). 

❖ Effect of time (t) : The velocity profile for different values of time ( t = 0.2, 0.4, 0.6, 0.8 ) is shown in 
Figure 6. It is clear that the velocity profile decreases with increasing time ( t ). 

❖ Effect of Eckert Number ( Ec ) : The Temperature profile for different values of Eckert Number 

(Ec = 0.5, 2, 4) is represented in Figure 7 and it is shown that the Temperature profile 

increases with increasing Eckert Number ( Ec ). 

♦♦♦ Effect of prandtl Number ( Pr ) : The Temperature profile for different values of Prandtl Number 
( Pr = 0.71, 7.0, 11.4 ) is shown in Figure 8 and it is observed that the Temperature decreases with 
increasing Prandtl Number ( Pr ). 

❖ Effect of Schmidt Number ( Sc ) : The Concentration profile for different values of Schmidt Number 

(Sc = 0.22, 0.6, 0.78, 0.96 ) is presented in Figure 9. It is observed from the figure that the 

concentration decreases with increasing Schmidt Number ( Sc ). 
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V. CONCLUSIONS 

The effect of different parameters like, phase angle(cot), Thermal Grashof Number(Gr), modified 
Grashof Number(Gc), Permeability parameter(k), Eckert Number(Ec), Prandlt Number(Pr), Schmidt 
Number(Sc), and time(t) are shown through graphs. 

> It is clearly observed the the velocity profile increases when Thermal Grashof Number (Gr) and modified 
Grashof Number (Gc) increases. But velocity profile decreases when Phase angle (cot), Permeability 
parameter (k), time (t) increases. 

> It is clearly observed the Temperature profile increases when Eckert Number (Ec) increases. But 
Temperature profile decreases when Prandlt Number (Pr) increases. 

> It is clearly observed the concentration profile decreases when Schmidt Number (Sc) increases. 




Fig 2: Velocity profile for different values of ( Ot ) Fig 4: Velocity profile for different values of GC 
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Fig- 6 : Velocity profile for different values of t Fig 3: T eiuperatiire profile for different values of Pr 
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Fig 7: Temperature profile for different values [if Ec 
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